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Abstract. The percolation phase transitions of two-dimensional lattice networks under a generalized 
Achlioptas process (GAP) are investigated. During the GAP, two edges are chosen randomly from the 
lattice and the edge with minimum product of the two connecting cluster sizes is taken as the next occu- 
pied bond with a probability p. At p = 0.5, the GAP becomes the random growth model and leads to the 
minority product rule at p = 1. Using the finite-size scaling analysis, we find that the percolation phase 
transitions of these systems with 0.5 < p < 1 are always continuous and their critical exponents depend 
on p. Therefore, the universality class of the critical phenomena in two-dimensional lattice networks under 
the GAP is related to the probability parameter p in addition. 


PACS. 64.60.ah 64.60.De — 89.75.Da 89,75.Hce 


1 Introduction 


The percolation phase transition concerns the formation of 
a macroscopic component in systems on both lattices and 
networks [i]. It provides a model for the onset of a macro- 
scopic component in random media [I] and social networks 
[2]. It was widely believed that the percolation transition 
is a typical continuous phase transition for various net- 
works [3]. However, Achlioptas, D’Souza, and Spencer [4] 
found recently that the percolation phase transition in 
random network becomes discontinuous (first-order) un- 
der the Achlioptas process (AP), where the edge with min- 
imum product of cluster masses is connected from the two 
randomly chosen unoccupied edges. The Achlioptas pro- 
cess suppresses the appearance of larger cluster and dis- 
courages the formation of a giant component, which has 
the size comparable with the number of vertices N. The 
percolation phase transition is delayed by this Achlioptas 
process and becomes sharper. It was argued by them|4] 
that the phase transition is discontinuous and named as an 
explosive percolation. Later, the Achlioptas process was 
introduced to the two-dimensional regular lattice networks 
and scale-free networks [6/7]. It was claimed that the 
explosive percolation was found both in lattice [5] and in 
scale-free networks [6][7]. 

In the article of Achlioptas et al.{4], the step interval 
A between the size of the largest component Sı = N12 
and Sı = 0.5N is used as the criterion of continuous or 
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discontinuous phase transition. Later, Ziff applied this cri- 
terion to the two-dimensional regular lattice networks|[5] 
[8]. It was found that the size-dependence of A in lattice 
network is quite different from that in random network. 
It is not well established that the first-order phase transi- 
tion can be distinguished from the continuous phase tran- 
sition by the size-dependence of A. It was argued by da 
Costa et al.[9] that the explosive percolation transition, 
under their modified Achlioptas process, is actually con- 
tinuous. Recently, Riordan et al.{10] show mathematically 
that all Achlioptas process have continuous phase transi- 
tions. The finite-size behavior of the order-parameter dis- 
tribution function has been used as the evidence of both 
discontinuous[L]] and continuous[12] phase transition. It 
is also argued with the finite-size scaling that the explosive 
percolation is continuous [[3\74]. This controversy about 
the character of explosive percolation is going on and calls 
for further investigations. 


In this paper, we investigate the percolation phase 
transition in two-dimensional lattice network under a gen- 
eralized Achlioptas process (GAP), which will be intro- 
duced in the next section. The generalized Achlioptas pro- 
cess is characterized by a probability parameter p. The 
GAP becomes the random growth model at p = 1/2 and 
the minority product rule at p = 1. Using the finite- 
size scaling analysis, our Monte Carlo simulation results 
demonstrate clearly that the percolation phase transition 
in two-dimensional lattice networks under the GAP is con- 
tinuous. It will be shown that the critical exponents and 
therefore the universality class of the continuous percola- 


2 Maoxin Liu et al.: Title Suppressed Due to Excessive Length 


tion phase transition depend on the probability parameter 
D. 

Our paper is organized as follows. In the next sec- 
tion, we introduce a generalized Achlioptas process in two- 
dimensional lattice network. In Section 3, we investigate 
the critical points of two-dimensional lattice network un- 
der the GAP and their critical exponents with the use of 
finite-size scaling. In Section 4, the finite-size scaling func- 
tion of the ratio S2/S1 is obtained at different probability 
parameter p, where S2 and Sı are the size of the second 
largest and the largest cluster in the network. The univer- 
sality class of the critical points in our model is discussed 
in Section 5. Finally we make some conclusions in Section 
6. 


2 Two-dimensional lattice network under the 
GAP 


We consider a two-dimensional square lattice with size 
L x L and periodic boundary conditions in both direc- 
tions. There are N = L? vertices in this lattice. We in- 
troduce a generalized Achlioptas process for adding edges 
into this lattice. In the generalized Achlioptas process, two 
edges are picked up randomly at each step. Each edge is 
connected with two clusters. The edge with the minimum 
product of the cluster sizes is chosen and added into the 
lattice with a probability p, where 0 < p < 1. Correspond- 
ingly, another edge is chosen with a probability 1 — p. At 
p = 0.5, the GAP is equivalent to the classic Erdés-Rényi 
(ER) rule, where edges are picked up randomly. The two- 
dimensional square lattice with the ER rule is actually the 
two-dimensional bond percolation (BP) model. The GAP 
at p = 1 is the product rule (PR) of Ref.[4] and our model 
becomes the PR model on the two-dimensional regular 
lattice. 

In our Monte Carlo simulations, there are N isolated 
vertices in a two-dimensional lattice at the beginning and 
then edges are added into the lattice through the GAP. 
With the edges added, we obtain a network in the lattice. 
The lattice network can be characterized by a reduced 
edge number r = N,./N, where N, is the number of the 
edges added. 

We have used the algorithm of Newmann and Ziff 
[I6] in our Monte Carlo simulations. For the investigations 
related only to the largest cluster of lattice networks in 
Figs. 3, 5 and 7, the linear sizes L = 32, 64, 128, 256, 512, 
and 1024 are taken. When the second largest cluster in 
the lattice is taken into account in addition, only three 
linear sizes L = 64, 128, and 256 are taken in Figs. 2, 
4 and 6. To get enough samples for the average of each 
simulation, different steps are taken for the simulation of 
different system size. In our Monte carlo simulations, we 
run 10,000,000 steps for L = 32 until to 6, 400, 000 steps 
for L = 1024. 

For the cluster ranked R and with size Sp(r, L; p), we 
defined its reduced size as 


sr(r, L; p) = Sr(r, L; p)/N. (1) 


Fig. 1. Reduced size of the largest cluster, defined in Eq. (I), as 
a function of the reduced edge number r at the lattice size L = 
1024 and probability parameters p = 0.5, 0.6, 0.7, 0.8, 0.9, 
and 1.0. 


In Fig [I] we shown the reduced size s1(r, L; p) of the largest 
cluster. 

In the Monte Carlo simulations of these data, we take 
the lattice size L = 1024 and the probability parameter 
p = 0.5, 0.6, 0.7, 0.8, 0.9, and 1.0. At small r, the re- 
duced size of the largest cluster is nearly zero. When r is 
large enough, the reduced size s; becomes finite. This in- 
dicates the formation of a macroscopic component. There- 
fore, there is a percolation phase transition in the lattice 
network. The transition value of the reduced edge number 
Te depends on the probability parameter p. It is shown in 
Fig fi] that re increases with the probability parameter p. 
This is plausible since a larger p means stronger suppres- 
sion of larger cluster and therefore the later appearance 
of a macroscopic component. In the following, we will try 
to verify that these percolation phase transitions are con- 
tinuous or not. 


3 Critical points of two-dimensional lattice 
networks under GAP 


If the percolation phase transitions above were continuous, 
the reduced size of the cluster ranked R should follow the 
finite-size scaling form [17/18] 

sr(r, L; p) = LP šR(tL”; p), 


(2) 


where t = (r — re)/re characterizes the deviation from the 
critical point re and v is the critical exponent of the cor- 
relation length € = ĉoļt| 7”. This finite-size scaling form 
is supposed to be valid in the asymptotic critical region 
where L > lattice spacing and |t| « 1. Outside the asymp- 
totic region, additional correction terms should be taken 
into account. 

For the largest cluster of lattice network, we have the 
finite-size scaling form 


sı (r, L; p) = LPs (tL; p). (3) 
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In the bulk limit L — ov, the reduced size of the largest 
cluster becomes 
si(r, 00; p) =0 (4) 


for r < re and 
s1 (r, 00; p) = ap t’ (5) 


for r > re. The emergent macroscopic component is char- 
acterized by the critical exponent 8. The smaller the crit- 
ical exponent ( is, the larger is the macroscopic compo- 
nent. 

Near the critical point, the reduced size of the second 
largest cluster can be written also in a finite-size scaling 


form 
s2(r, L; p) = LP ša (tL; p). (6) 


Using Eqs. (3), and (6), we can obtain the finite-size 
scaling form of the ratio 


S2/S1 = 5o(tL/”; p)/51 (tL; p) = UL”; p). (7) 
At the critical point t = 0, the ratio 
52/Si|,-,, = U (0; p), (8) 


which is independent of the system size L. Therefore, the 
curves of 52/5; at different system size L have a cross- 
point at r = re. The critical point corresponds the fixed 
point of S2/S1, which can be used to determine the critical 
point of our system. 

The logarithm of Eq. (8) can be expressed as 


ln sı(r, L; p) = —(8/v)n L +n (tL”; p). (9) 
At the critical point r = re, we have 
In sı(re, L; p) = —(8/v) ln L + Ins} (0; p), (10) 


which is a straight line with respect to In L. We can use 
this property to determine the critical point re of the lat- 
tice networks also. From the slope of this straight line, the 
critical exponent ratio 8/v can be determined. 

In the following, we use both the fixed point of S2/91 
and the linear dependence of In sı on In L as the criterion 
to determine the critical point of the two-dimensional lat- 
tice networks under the generalized Achlioptas process. If 
we have reached the asymptotic critical region, both sı 
and sz satisfy the finite-size scaling form in Egqs.(3) and 
(6). The critical reduced edge numbers re obtained from 
S/S and Ins; should be equal. On the other hand, the 
consistence of re obtained from two different methods can 
be used as the indicator of the accuracy of our simulation 
results. 

For the generalized Achlioptas process with probabil- 
ity parameter p = 0.5, edges are added randomly into a 
two-dimensional lattice and our model becomes the so- 
called bond percolation model. It is well known that the 
bond percolation model has a continuous phase transi- 
tion. We can determine the critical point of this model 
by the two methods described above. In Fig] the ratio 
S/S; is shown as a function of the reduced edge num- 
ber r at different system size L. A fixed point between 
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Fig. 2. Size ratio 52/51 of the second largest to the largest 
cluster at p = 0.5 and different L. There is a fixed point at 
Te = 0.5000 + 0.0002. 
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Fig. 3. Log-log plot of the reduced size sı of the largest cluster 


at p = 0.5 and different r. We have chosen linear sizes L = 
32, 64, 128, 256, 512, and 1024. 


r = 0.4998 and r = 0.5002 is found and is in full agree- 
ment with re = 1/2 of the bond percolation model. In 
Table[i] we denote the critical point obtained from S2/S1 


as rí} = 0.5000 + 0.0002. 


As we have discussed above, a critical point can be de- 
termined alternatively by the linear relationship between 
ln sı and ln L. In Fig] In sı (r, L; p) at different reduced 
edge numbers are shown. At r = 0.4996 and r = 0.5004, 
the curves of In sı (r, L; p) are curved and their curvatures 
have different sign. At r = 0.5, the curve of ln sı (r, L; p) 
becomes a straight line. Therefore, we obtain the criti- 
cal point r? = 0.5000 + 0.0004, which is in agreement 
with r®. From the slope of ln sı(r, L; p) at re, we get 
the critical exponent ratio 8/v = 0.108. Our Monte Carlo 
simulation results agree very well with the exact results 
re = 1/2 and 8/v = 5/48 of two-dimensional bond perco- 
lation model, for a review, see Ref. [I9]. 
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Fig. 4. Ratio S2/S; at p = 0.8 and different L. There is a 
fixed point at r = 0.5208 + 0.0001. 


For probability parameter p = 0.8, the edge that mini- 
mizes the product of two connecting cluster sizes is added 
into the lattice with a probability 0.8 from two randomly 
chosen edges. The connection of smaller clusters is fa- 
vored. The critical point of this system is investigated 
by the fixed point of S2/Sı and the linear dependence 
of Ins; on ln L. In Fig] the ratio S2/51 is plotted as a 
function of the reduced edge number r at different sys- 
tem sizes L. A fixed point of 52/51 is found. It is be- 
tween r = 0.5207 and r = 0.5209. So there is a continuous 
phase transition in this system and the critical point is at 
ph) = 0.5208+0.0001. Alternatively, this critical point can 
be determined from In sı (r, L; p). In Fig] In s1(r, L; p) at 
r = 0.5205, 0.5207, and 0.5209 are shown. The curvature 
of Insi(r, L; p) is negative at r = 0.5205 and positive at 
r = 0.5209. The function In sı (r, L; p) at r = 0.5207 can be 
described quite well by a straight line with zero curvature. 
So the critical point rP? = 0.5207 + 0.0002, in agreement 
with ri), The slope of the straight line at re = 0.5207 
gives the critical exponent ratio 8/v = 0.081. 

At p = 1.0, our model becomes the two-dimensional 
lattice network under the Achlioptas process. The ratio 
S/S, of this model is shown in Fig[6] for different sys- 
tem size L. There is a cross-point between r = 0.5265 and 
r = 0.5267, which corresponds to the critical point of this 
system. Therefore its critical point is at r® = 0.5266 + 
0.0001. The curves of In sı (r, L; p) at r = 0.52651, 0.52655, 
and 0.52659 are shown in Fig[7] Their curvatures change 
with r from negative to positive. The function becomes 
linear with respect In L around r = 0.52655. So we get the 
critical reduced edge number r? = 0.52655 + 0.00004, 


which agrees with r{® given above. From the slope of 
Insi(r-, L; p) with respect to ln L, the critical exponent 
ratio B/v = 0.064 is obtained. Our results of the critical 
point and the critical exponent ratio are in full agreement 
with the results of Refs. [51617]. 


In Fig [8] we summarize the critical reduced edge num- 
bers re of two-dimensional lattice networks under the GAP 
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Fig. 5. Log-log plot of the reduced size sı of the largest cluster 
at p = 0.8 and different r. We have chosen linear sizes L = 
32, 64, 128, 256, 512, and 1024. 
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Fig. 6. Ratio S2/S1 at p = 1.0 and different L . There is a 
fixed point at r = 0.5266 + 0.0001. 
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Fig. 7. Log-log plot of the reduced size sı of the largest at 
p = 1.0 and different r. We have chosen linear sizes L = 
32, 64, 128, 256, 512, and 1024. 
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Fig. 8. Dependence of the critical reduced edge number re on 
the probability parameter p. The error bars of Monte Carlo 
simulation data are smaller than the symbol. 
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Fig. 9. Dependence of the critical exponent ratio B/v on the 
probability parameter p. 


with different p. It is found that re increases with p. At a 
larger p, the continuous percolation phase transition ap- 
pears at a larger critical reduced edge number r, and the 
formation of a giant component is delayed. 

In Fig[9] the dependence of the critical exponent ra- 
tio 8/v on p is shown. With the increase of p, the ratio 
G/v decreases. We will show in the next section that 1/v 
increases with p. So it can be concluded that the criti- 
cal exponent 8 decreases with p. Smaller 8 indicates the 
stronger emergence of a giant component in the networks 
after the percolation transition. Therefore, the increase of 
p results in the delayed appearance of a continuous perco- 
lation phase transition and the formation of a larger giant 
component at the same time. 


4 Finite-size scaling functions of S/S; 


In the last section, we have mentioned that the size ratio 
S2/Sı follows the finite-size scaling form in Eq. (Z) when r 
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Fig. 10. Finite-size scaling function U(tL‘/”; p) of the ratio 
S2/Si at p = 0.5. The critical exponent of two-dimensional 
bond percolation v = 4/3 is taken for the scaling variable. 


is near the critical point re. For a given probability param- 
eter p, the different curves of S2/S; at different L collapse 
into a finite-size scaling function after using the scaling 
variable tL!/”, where v is the critical exponent of cor- 
relation length. In the following, we will investigate the 
finite-size scaling function of S2/S1 for different p. 


At p = 0.5, we use the critical exponent of two-dimensional 


bond percolation model v = 4/3 [I9] for the finite-size 
scaling function of $2/S,. After defining the scaling vari- 
able tL!/” with this value of v, our Monte Carlo simulation 
results at L = 64, 126, 256 collapse and we get the finite- 
size scaling function of S2/S,, which is shown in Figo] 

At p = 1.0, the critical exponent v is unknown. Ac- 
cording to the finite-size scaling form in Eq. (Z), the curves 
of S/S) at different system sizes can collapse only when 
the correct critical exponent v is used for the scaling vari- 
able. This property can be used also for determining the 
critical exponent v. At 1/v = 0.93, the curves of S/S} 
at L = 64, 128, 256 collapse into its finite-size scaling 
function, which is shown in Fig [ii] 

In Fig[12] we demonstrate the variation of the finite- 
size scaling function of S2/S1 with the probability parame- 
ter p. In the region before the percolation phase transition, 
the finite-size scaling function of 2/5) increases with p. 
The second largest cluster in this region is more important 
at larger p. In the region after the percolation phase tran- 
sition, the finite-size scaling function of S2/S, decreases 
with p. The largest cluster in this region is more domi- 
nant at larger p. To get the finite-size scaling functions 
at different p, the corresponding exponents of correlation 
length are determined and presented in Fig [13] and Tabel 


5 Universality classes 


The concept of universality plays a fundamental role in 
statistical and elementary particle physics [20)[2i]. The 
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Fig. 11. Finite-size scaling function U(tL™”; p) of the ratio 
S2/S1 at p = 1.0. The inverse of the critical exponent 1/v = 
0.93 is taken for the scaling variable. 
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Fig. 12. Finite-size scaling functions U(tL'/”; p) of the ratio 
S2/S1 at p = 0.5, 0.6, 0.7, 0.8, 0.9, and 1.0. The corresponding 
1/v for the scaling variable at different p are presented in Fig[I3] 
and Tabel [I] 


Table 1. Critical reduced edge number re, ration of critical 
exponent 3/v and inverse of critical exponent 1/v at different 
probability parameters. We obtain the critical reduced edge 
number r® from S2/Sı and r2 from ln sı. 


p rp rP B/v 1/v 

0.5 0.5000(2) 050004) 01084) 0.75 

0.6 0.5082(2) 0.5082(4) 0.102(3) 0.77(1) 
0.7 0.5153(3) 0.5153(2) 0.092(4) 0.79(1) 
0.8 0.5208(1) 0.5207(2) 0.081(5) 0.83(1) 
0.9 0.5244(1) 0.5244(1) 0.070(3) 0.88(1) 
1.0 0.5266(1) 0.52655(4) 0.064(3) 0.93(1) 
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1/v 


Fig. 13. p-dependence of the inverse of critical exponent 1/v. 
The Monte Carlo data at p = 0.5, 0.6, 0.7, 0.8, 0.9, and 1.0 
are shown. 


Fig. 14. p-dependence of the finite-size scaling function 
U(tL/”; p) of the ratio S2/S; at the critical point t = 0. 


universality is characterized by the dimensionality d of 
the system and by the number n of the components of 
the order parameter [22]. Within a certain (d, n) univer- 
sality class, the critical exponents are independent of mi- 
croscopic details and are universal. In a finite-size system 
near its critical point, there is also universality. For ex- 
ample, the Binder cumulant ratio of magnetization at the 
critical point is universal. The ratio S2/S1 here is similar 
to the Binder cumulant ratio of magnetization. We could 
suppose that the ratio $2/S; at the critical point is also 
universal and does not depend on microscopic details. 

In our previous investigations of the two-dimensional 
lattice networks under a generalized Achlioptas process, 
the dimensionality d of systems is fixed and the macro- 
scopic property of their order parameter is unchanged. 
But we have found in Figs[9[13] and [4] that the critical 
exponents 8, v and the ratio S2/S; at re depend on the 
probability parameter p. So the universality of percola- 
tion phase transition in these networks is characterized in 
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addition by the probability parameter p of GAP. A dif- 
ferent probability parameter p of GAP generates a differ- 
ent probability distribution of configuration. So the prob- 
ability parameter p is actually related to the macroscopic 
character of network. For a general classification of uni- 
versality class in complex networks, further investigations 
are needed. 


6 Conclusions 


We have investigated the percolation phase transitions 
in two-dimensional lattice network under a generalized 
Achlioptas process. In this GAP, we choose randomly two 
unoccupied edges in a two-dimensional lattice and the 
edge that minimizes the product of the two connecting 
cluster sizes is taken with a probability p. Our model 
becomes the two-dimensional bond percolation model at 
p = 0.5 and the two-dimensional lattice network under the 
minority product rule at p = 1. 

The size Sı of the largest cluster in the lattice increases 
with the edge number N,. When the reduced edge number 
r = N,/N is larger than a certain value re, Sı becomes 
comparable with the lattice size N = L?. At re, a gi- 
ant component emerges and there is a percolation phase 
transition. From the finite-size scaling analysis of Sı and 
the ratio S2/S1, we can conclude that this percolation 
phase transition is continuous at probability parameter 
0.5 <p < 1. The critical exponent ratio G/v can be deter- 
mined from the power-law behavior of Sı at re. To obtain 
the finite-size scaling function of the ratio S2/S1 from the 
Monte carlo simulation data of different L with the scal- 
ing variable tL!/” , the critical exponent of correlation 
length v can be fixed. We find that the critical reduced 
edge number re increases with the probability parameter 
p, which is shown in Fig[8] The critical exponent ratio B/v 
and the critical exponent v decrease with the probability 
parameter p, as demonstrated in Fig[] and [3] Under the 
GAP with 0.5 < p < 1, the formation of larger cluster is 
suppressed and this suppression increases with p. So the 
formation of a giant component should be delayed at a 
larger probability parameter p. This delay is accompanied 
then by the stronger emergence of the giant component, 
which is characterized by smaller £. It is plausible that re 
increases and decreases with p. The finite-size scaling 
functions of the ratio S/S) are given for different p in 
Figi 

Within a certain universality class characterized by 
the dimensionality of the system and by the number of 
components of the order parameter, the universal quan- 
tities (critical exponents, amplitude ratios, and scaling 
functions) of different systems are identical. For the two- 
dimensional lattice networks under a GAP we discuss here, 
the critical exponents 8, v and the ratio 2/5} at the criti- 
cal point depend on the the probability parameter p, which 
has been pointed out above. So the universality class of 
the percolation phase transition in this model should be 
characterized in addition by the probability parameter p. 
To understand the universality class of the critical phe- 
nomena in networks in general, further investigations are 


needed. For random networks, we introduce also a gen- 
eralized Achlioptas process and the phase transitions in 
these systems are investigated [23]. 

This work is supported by the National Natural Sci- 
ence Foundation of China under grant 10835005. 
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